Nonlinear Kinetic Development of the Weibel Instability 
and the generation of electrostatic coherent structures 
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The nonlinear evolution of the Weibel instability driven by the anisotropy of the electron distri- 
bution function in a collisionless plasma is investigated in a spatially one-dimensional configuration 
with a Vlasov code in a two-dimensional velocity space. It is found that the electromagnetic fields 
generated by this instability cause a strong deformation of the electron distribution function in 
phase space, corresponding to highly filamented magnetic vortices. Eventually, these deformations 
lead to the generation of short wavelength Langmuir modes that form highly localized electrostatic 
. structures corresponding to jumps of the electrostatic potential. 
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I. INTRODUCTION 



The generation of macroscopic magnetic fields in a plasma is a fundamental process that occurs in laboratory and 
astrophysical plasma configurations in a wide range of regimes of plasma collisionality. This process appears naturally 
in a plasma. In plasma dynamics the "mechanical" and "electrodynamic" degrees of freedom couple very efficiently 
allowing for a rapid transfer from mechanical to electrodynamic energy and vice versa, as exemplified by the the 
dynamo mechanism [l| which is responsible for the generation of the magnetic fields of the planets and of the stars. 



Here we are interested in investigating the different physical mechanisms that can accompany the generation of 
. ■ magnetic fields in a collisionless plasma regime, i.e., under conditions that are very far from thermodynamic equi- 
librium, and in particular in the generation of longitudinal electric fields and of electrostatic potential jumps due to 
^ nonlinear effects and to the related modification of the electron distribution function in phase space that accompany 

p> ' the generation of the magnetic field. 

!— i. We refer in particular to the basic mechanism identified more than half a century ago by Weibel Q where the 
magnetic field is generated by the development of an instability driven by the temperature anisotropy of the electron 
distribution function in a collisionless plasma. Later on, a systematic linear kinetic analysis of the instabilities at 
play in an anisotropic plasma, without an initial mean magnetic field, has been carried out in an electron-ion plasma 
|3|]. We recall that other mechanisms of magnetic field generation can be at play in different collisionality regimes, as 
(f") . extensively reviewed in Ref. 

A rough physical model of the magnetic field generation can be formulated in terms of a "thermal" machine working 
qq , between the two 'temperatures', T x and T y which correspond to the two 'temperatures' that characterize the electron 
distribution function: T y describes the distribution width in v y and T x < T y the distribution width in v x . In a highly 
collisional regime the two temperatures would be rapidly equalized by the effect of collisions. In a collisionless plasma 
regime work can be extracted from this 'temperature' difference. That this work is associated with the generation of 
magnetic field, i.e. in other words that this work corresponds to the transformation of thermal energy into magnetic 
energy, can be guessed by means of a virtual displacement argument borrowed from the theory of the closely related 
current filamentation instability We imagine splitting at each point the electron distribution function into two 
parts, corresponding to positive and to negative values of v y respectively, and displacing them in the transverse plane 
with a virtual displacement of the form sin (k x x) = —£q x sin (k x x), where the upper indices refer to the populations 
with positive and with negative value of v y , respectively. In this way opposite current densities along y are formed 
in the plasma. These currents are modulated along x and produce a magnetic field along z. Since opposite currents 
repel, the initial displacement is reinforced, the instability can develop and the magnetic field along z can grow. 

Actually, the nonlinear evolution of the Weibel instability is far richer, even when we focus our investigation on 
a simplified geometrical configuration where all quantities are assumed to depend only on x and velocity space is 
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reduced to the v x -v y plane. In this context we recall that the restriction to such a 1D-2V configuration excludes, 
among other phenomena that can accompany the nonlinear evolution of the Weibel instability, the development of a 
secondary magnetic field line reconnection instability of the type discussed in Ref. @ . 

In the one dimensional configuration considered here, the development of the Weibel instability is accompanied by 
the generation of electric fields: an inductive electric field in the y direction intrinsically related to the variation in 
time of B z and a longitudinal field along x. For the sake of clarity, we recall that in the case of the related current fila- 
mentation instability, where an effective anisotropy in velocity space arises from the presence of two counterstreaming 
electron beams [5|, the growth of a longitudinal electric field in a two dimensional spatial configuration is due to the 
onset of a primary two-stream instability with wave vector along the direction of the beams (i.e. along y with the 
choice of axes considered here) . On the contrary, in the present case, the growth of a longitudinal electric field can 
only occur as a secondary process produced by the development of the Weibel instability itself. We also note that 
the generation of a longitudinal field, orthogonal to the beam direction, during the nonlinear evolution of the current 
filamentation instability in a one dimensional configuration was discussed analytically^, H| in terms of a two fluid, 
cold plasma description and shown to occur near wave breaking conditions due to the nonlinear coupling between the 
current filamentation instability and Langmuir waves. 

In the present article we recover within a kinetic description the nonlinear generation of a longitudinal field, 
orthogonal to the direction of the perturbed current produced by the Weibel instability in its linear phase. Furthermore 
we show that the nonlinear time evolution of the electron distribution function in phase space, in the presence of the 
magnetic and inductive electric fields produced by the Weibel instability, is characterized by a mixing process that 
occurs mostly at the position along x where the amplitude of the magnetic field is largest and by particle acceleration 
or deceleration that occurs mostly at the position along x where the amplitude of the inductive electric field is largest. 
The combination of these processes leads to the formation of a highly non monotonic distribution function in velocity 
space that can further excite Langmuir waves resonantly. Eventually, the evolution of the longitudinal fields is found 
to lead to the formation of multipolar electrostatic structures characterized by a bi-polar or tri-polar electric field 
profile that are reminiscent of the "isolated electrostatic structures" detected by satellites in- situ in the solar wind 
plasma, in the auroral regions and in the Earth magnetosphere and bow shock. These structures have a typical width 
of the order of a few tenth of a Debye length and are thought to propagate along the magnetic field lines at speeds 
smaller than the solar wind speed. The tri-polar structures are also known as double-layers and are accompanied 
by a net potential drop. They have recently been proposed as a possible candidate in order to account, through a 
sequence of such double layers, for the total potential difference from the solar corona to the Earth, as consistent with 
exospheric models (see the recent review [9j and references therein). 



II. GOVERNING EQUATIONS 

We consider high frequency modes evolving on electron time scales defined by the inverse of plasma frequency 
in a collisionless plasma and describe the dynamics of the plasma using the Vlasov-Maxwell system of equations for 
the electron distribution function f e taking the ions to remain at rest. Thus in dimensionless units we write 

f +v.^ + ( E + vxB). f=0, (1) 
at ox ov 

where velocities are normalized to the speed of light and times to the inverse of the plasma frequency uj pe . Lengths 
are thus normalized to the electron skin depth d e = c/uj pe . 

We restrict ourselves to a 1D-2V configuration where all quantities depend on x and time only, the particle velocities 
and the electric field have x-y components and the magnetic field is along the z axis, and write Maxwell's equations 
in the form 



dB z dE y dB z _dE y , d 



2, 



J v TaZ = ~P> ( 2 ) 



dt dx ' dx dt d 2 x 

where p and j are the dimensionless charge and current densities and 4> the dimensionless electrostatic potential. 

The Vlasov equation (TTJ) is integrated in the phase space (x, v x , v v ) with periodic boundary conditions along the x 
direction. For a description of the adopted code see Ref. We take the following initial condition 

N 

fe(x,v x ,v y ,t = 0) = f M (v x ,v y )[l + cos(k n x + <p n ))}; k„ = 2nn/L x . (3) 
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FIG. 1: Left frame: amplitude of the Fourier components of the magnetic field B Zt k at t = 300 vs. k. Right frame: time 
evolution of the Fourier components of the most unstable k — k max = 1 magnetic and inductive electric field mode, B z ^=i and 
E V: k=i, solid and dashed line respectively, and of the most unstable k = 2.3 longitudinal electric field E x ,k=2.3 (dotted line). 
The dash-three dotted line represents the growth rate of the magnetic field. 



Here L x is the length of the simulation domain in the x direction, e is the perturbation amplitude varied in the range 
10~ 6 < e < 10~ 4 and Jm(v) is a two-temperature Maxwellian distribution given by 

71 / iP 1 V \ 

fu(v x ,v y )= J-^ expi-^-^), (4) 

with T x = 4 • 10 -4 corresponding to a plasma with a non-relativistic temperature, T y /T x = 12 and </>„ are random 
phases. At t = we also introduce a perturbation on the magnetic field: 

N 

B z = a%y^ J [cos(k n x + i/j n )}, (5) 

71=1 

where is the initial field amplitude for each wave and ip n are random phases. We take = 1CP 4 and N = 200 
in both Eqs. ((3]) and ([5]). We note the initial conditions with random initial amplitudes of the density and magnetic 
and electric field have been also used without any significative qualitative and quantitative change in the results. 

We consider a field free plasma, hence < B z >= at t = 0. The computational velocity phase space is given 
by [— v™ ax , v™ ax ] x [— u™ oa: , v™ ax ]. We take v" lax — 0.1 and v " lax = 0.27 since the system is characterized by 
an initial temperature anisotropy T y > T x . However, the distribution function will spontaneously evolve 
towards a more isotropic configuration by reducing its width in the v y direction and by generating small 
scale filamentary structures with the same characteristic size in both directions. Therefore, we take a 
larger number of points in the v y direction in order to reduce the difference between the resolution in 
the two directions. For computational reasons (i.e. not to increase the total number of grid points too 
much), a good compromise is N Vx = 281 and N Vy = 121. The time step is At = O.OOSa;^ 1 . The evolution of 
the system is investigated up to t = 500OC, 1 . The length of the spatial simulation domain is L x = 6ird e . 
The space discretization along x is dx = 0.03 d e , where d e is the electron skin depth, corresponding to 
N x — 600 grid points in the x direction. This allows us to account for phenomena that occur on small 
spatial scales that are larger than, but of the order of, the electron Debye scale Xd = 0.02 d e . 



III. ONSET OF THE WEIBEL INSTABILITY AND EARLY NONLINEAR PHASE 



The onset and the linear growth of the Weibel instability up to the beginning of its saturation phase is shown in 
Fig. [TJ as obtained from the numerical integration of the Vlasov equation described above. In the left frame the 
Fourier amplitude of the magnetic field B z ^ at t = 300 is shown vs. k (note that only a part of the interval in fc is 
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FIG. 2: Growth of k — 1,2,2.3 Fourier components of the longitudinal field E Xt t, dashed, dotted and solid line, respectively. 

shown). We see that the modes in the range 0.3 < k < 1.5 have grown significantly with their maximum growth rate 
corresponding to k = 1, i.e. in dimensional units to k — d~ x . In the right frame, the amplitude of the magnetic 
field B Zj k vs. time is plotted for the most unstable mode, k = 1 (solid line). After the initial transient phase the 
amplitude of the magnetic held grows exponentially, for t < 270 (linear phase), with growth rate ~/ m ax ~ 0.04 (see 
dash-three dotted line), consistent with the analytical results (see Ref. Q). 

The growth of the inductive electric field component E y j, associated to B z ^ is shown in the same frame (dashed 
line) still for k = 1 , together with the longitudinal field E x ^ (dotted line) taken at k = 2.3. We note the E Vt k 
reaches its maximum value at the saturation time of B z ^ and then starts to decrease. This is consistent with the fact 
that the ratio between the inductive electric field and the magnetic field generated by the Weibel instability scales 
with the instability growth rate so that, as soon as the Weibel instability saturates, the mode structure is dominated 
by the magnetic field. In this phase the longitudinal electric field E Xi j~ arises because of the coupling between the 
Weibel instability and the Langmuir waves, as discussed within the fluid approximation in Ref. Q, due to the electron 
density modulation induced by the spatial modulation of B^ k , and thus grows in time at twice the growth rate of the 
magnetic field. 

While the Fourier component of the longitudinal field shown in Fig. [T] right frame, corresponds to k = 2.3, the 
spectrum in k of the longitudinal field and of the density modulation is wider, as consistent with the fact that the 
Weibel instability grows with very close growth rates in a relatively large range of values of k around k ~ 1. As shown 
in Fig. [2J different Fourier components of the longitudinal field grow with nearly equal values of the growth rate. As 
a consequence, while the longitudinal field grows two times as fast as the magnetic field, its spatial structure contains 
additional wavelengths besides A/2 = it. 

The frequency spectra of B z ^ for fixed k, are shown in Fig. [3] for k = 1,8, 10. The spectral amplitude is largest 
for k = 1 and in the low frequency part of the frequency spectrum corresponding to the Weibel instability. For 
larger values of k the frequency spectra exhibit a weaker feature for large frequencies corresponding to low amplitude 
transverse electromagnetic waves propagating essentially at the speed of light (to ~ k in dimensionless units) . 

The frequency spectra of the longitudinal field E x ^ at fixed k, are shown in Fig. 0]for k = 1,2,8 and exhibit a 
clear peak at the electron plasma frequency. This peak appears to be less pronounced for larger values of k. 

Additional numerical integrations of the Vlasov equation, not shown here, for different initial 
anisotropy ratios (T y /T x = 9 and T y /T x = 4) at constant T x indicate that, while the k = 1 mode re- 
mains the most unstable electromagnetic mode, the spectrum of the excited modes becomes narrower 
with decreasing anisotropy, since higher k modes grow less fast, as consistent with linear theory. This 
effect is also present in the spectrum of the longitudinal electric field E x j. in which case the value of k 
corresponding to the maximum growth rate decreases as the anisotropy ratio is decreased. The maxi- 
mum value of B z , taken at the time when its growth is saturated, decreases with decreasing anisotropy 
ratio. Conversely, if T x is decreased at constant anisotropy ratio (T y /T x = 12), the maximum value 
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FIG. 3: Frequency spectra of the magnetic field component B z> k for k — 1, 8, 10 (from left to right). 
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FIG. 4: Frequency spectra of the longitudinal electric field component E Xl k for k = 1, 2, 8 (from left to right). 



of B z is almost unchanged, while the value of k corresponding to the maximum growth rate of the 
longitudinal electric field decreases from k = 2.3 at T x = 1 to k = 2 at T x = 0.5. 



IV. PHASE SPACE EVOLUTION IN THE ADVANCED NONLINEAR PHASE 



During the linear growth of the instability, the evolution of the electron distribution function in velocity space is 
characterized by a differential rotation in velocity space at the points where the absolute value \B Z \ of the magnetic 
field generated by the Weibel instability has a maximum and by a Y- shaped deformation with axis along v y at the 
points where the perturbed magnetic field vanishes which correspond to the points where the absolute value \E y \ of 
the inductive electric field is largest. This behaviour is shown by the frames in the second row of Fig. [5] The rotation 
of the distribution function is clockwise or anti-clockwise depending on the sign of B z (labelled as B max and B min 
respectively in Fig. [5j while the Y-shaped deformation is turned upside down depending on the sign of E y (only one 
sign is shown in Fig. [5]). As the Weibel instability enters its fully nonlinear phase the winding of the distribution 
function becomes tighter and the Y- deformation more marked until both features become " multi-armed" , as shown by 
the last row in Fig. [5] The winding of the distribution function and the formation of phase space vortices associated 
to the magnetic field B z {x) is also apparent from the contour plots of the distribution function at t — 340 in the x-v x 
space shown in Fig. [5] for different values of v y . 

Disregarding for the moment the effect of the longitudinal electric field, this complex behaviour can be attributed to 
the combined effects of the magnetic field B z , that deflects the electron orbits (where \B Z \ is largest and E y vanishes) 
and of the inductive electric field E y that accelerates or decelerates the electrons along y (where E y is largest and 
\B Z \ vanishes). In order to visualize this combined effect of the instability fields on the plasma electrons a set of 
particle orbits have been computed by direct integration of the particle equation of motion in the electromagnetic 
fields obtained from the Vlasov integration (shown in Fig. [7] in {x,t) space). 

For illustration a representative orbit is shown in Fig. [8] The particle is initialized at x = 10 with vox = 0.01 
and vq v — 0.05 and moves straight until around t — 270 the electromagnetic fields of the instability have grown 
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FIG. 5: Time evolution of the electron distribution function with t = 0(first row), t — 250(second row), t = 
290(third row) and t = 330(last row) at different positions along x = 7.8,6.4, 10.2 (from left to right). Red (blue) 
part is most (least) densely populated. 

significantly. At this time this particle undergoes a net deceleration and its kinetic energy decreases while vq x 
and vq v become oscillatory. The same oscillatory behaviour is displayed by the particle kinetic energy. The particle 
trajectory becomes trapped along x between two positions, separated by roughly half a wavelength of the dominant k 
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FIG. 6: Contour plots of the distribution function in the x-v x space at t = 340 for different values of v y 
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FIG. 7: Contour plots showing the combined x-t dependence of B z (left frame) and of E y (right frame). 

component of B z , where the amplitude of the magnetic field is maximum (B z is negative near x = 10.5 and positive 
near x = 12.5). 

As shown by the contour plots in Fig. [5] the electron distribution function becomes "multi-armed" as the instability 
evolves. This implies that the distribution function develops steep positive slopes for velocities of the order of the 
electron thermal velocity along x as shown in Fig. [9] Although these positive slopes evolve in time as the electron 
distribution function becomes increasingly twisted, they can give rise to a further excitation of Langmuir waves but 
with phase velocities much smaller than those of the Langmuir waves driven by the nonlinear coupling discussed in 
Sec. IIII1 This new destabilizing process is shown in Fig. [SJ last frame, where the time evolution of the longitudinal 
electric field component is shown versus time for k = 25 which corresponds to a phase velocity along x equal to 0.04. 
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FIG. 8: Example of an electron orbit in the time evolving electromagnetic fields: particle velocities vs. time (top frames), 
particle trajectory (bottom left frame), particle kinetic energy vs. time (bottom right frame). 

Contrary to the components of the longitudinal field with large phase velocities shown in Fig. [5J in the initial phase 
this high-fc, low phase velocity longitudinal field decays exponentially to noise level, as consistent with the strong 
Landau damping due to the initial Maxwellian shape of the electron distribution function. Later, around t = 280, the 
non monotonic feature of the electron distribution function shown in the first frame in Fig. [9] makes the mode grow 
resonantly with a growth rate of the order of 7 = 0.15, i.e. faster than the low-fc modes produced by nonlinear density 
modulations. Eventually the mode reaches an amplitude that is smaller by approximately an order of magnitude than 
that of the low k longitudinal modes. 

A similar growth is also exhibited by longitudinal modes within a wide range of values of k that extends towards 
larger values of k and thus towards smaller phase velocities that resonate well inside the electron distribution function. 
The growth of these lower phase velocity modes starts somewhat earlier than that of the k = 25 shown in Fig. [5] but 
reaches significantly lower final amplitudes. 

This destabilization mechanism shows that in a collisionless plasma dominated by kinetic effects the nonlinear 
energy transfer between modes with different wave-number can occur nonlocally in k space, with modes with higher 
values of k being excited earlier than modes with comparatively lower values of k. 

As mentioned before, smaller values of the magnetic field B z are found when smaller anisotropy ratios 
are considered. In this case the twisting of the electron distribution function is less pronounced and, 
consequently, the excitation of Langmuir waves with phase velocities resonating inside the electron 
distribution function is considerably weaker. 

V. FORMATION OF ELECTROSTATIC COHERENT STRUCTURES 

The wide spectrum of longitudinal modes that are excited by the phase space resonances discussed above leads to 
a highly structured electron density distribution n e along x, while the k spectrum magnetic field generated by the 
Weibel instability is considerably narrower and B z remains spatially regular even at late times (see Fig. [7J left 
frame). This is also illustrated in Fig. 1101 where we show the spatial profile of B z and E x versus x 
at t = 468. It is seen that the electric field, besides exhibiting a spatial modulation at roughly half the dominant 
magnetic field wavelength, has developed several narrow spikes, in particular a very narrow one at x ~ 13.6 
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FIG. 9: Electron distribution function at x = 7.8 and v y — vs. v x at t = 280, 310, 320 (first three frames). Time evolution of 
the longitudinal electric field component E X: k=2S (last frame). Electron distribution function at t — in the first frame with 
dashed lines. 




FIG. 10: Plot of £? z and E x , left and right frame versus x at t = 468. 



(corresponding to a sharp minimum in the the electron density). In the one-dimensional configuration 
considered here, these spikes correspond to "multi layer" electrostatic structures accompanied by 
a drop in the potential. A detailed analysis shows that these structures are not a true equilibrium 
between the electrostatic field and the pressure gradient, as is generally found for double layer coherent 
structures associated with phase space holes Ref. [11]. However, since these structures survive on time 
scales of the order of several ion periods, we still define them as coherent electrostatic structures. The 
time evolution of the longitudinal electric field is shown in the first four frames in Fig. [TT]vs. x and in the last frames 
in Fig. ITT1 where E x at x = 13.6 is plotted vs. t. A clear view of the formation of the electrostatic structures is given 
in Fig. [12] where we show the electrostatic field in the (x, t) plane. We see that after forming large scale structures on 
the magnetic (one-half) scale around t — 250, slowly propagating small scales electrostatic multipolar structures are 
formed on the Debye length scale. As they propagate, these structures interact with each other leading to coalescence 
or to the formation of new structures with similar characteristics. We observe that these structures do not 
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FIG. 11: Top four frames: Plot of E x versus iati= 471,473,477,480. Last two frames: Plot of E x versus t at 
x = 13.6 



form in the numerical simulations with smaller anisotropy ratios (T y /T x — 4). 

Electrostatic features of this type are encountered under space conditions and were first observed almost thirty years 
ago in auroral regions in the form of dipolar electric structures as well as double layers [Hj]. These "large amplitude 
solitary structures" are in general observed with the electric field aligned with the mean magnetic field direction 
(l3j . In the solar wind observations of dipolar structure have been reported starting from ten years ago in Ref. [l4| . 
From a theoretical point of view, there is no consensus about the origin of such structures. Accelerated particles, or 
beams, generated far away by some energetic events such as, for example, magnetic reconnection, are probably the 
best candidate for the energetic driver of the mechanism at the basis of the formation of multipolar structures [15l4l7| . 
In this context, it is interesting to note that plasma temperature anisotropies, a typical feature expected in many 
regions in solar wind plasmas, may represent a possible source for their formation through a multi stage process that 
involves a complex evolution of the electron distribution function under the influence of the electromagnetic fields 
generated by the Weibel instability. A study of plasma electron temperature anisotropy in the presence of a weak, 
mean magnetic field (see [l8[ for the linear analysis and [lj| for simulations) is in progress. 
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FIG. 12: Contour plots in the (x — t) plane of the electrostatic field E x . 
VI. CONCLUSIONS 

Pressure anisotropy is a common feature of laboratories and of space plasma and provides the free energy source for 
the development of the Weibel instability and the generation of a quasi static magnetic field. In this paper we have 
investigated the nonlinear evolution of the Weibel instability arising from the anisotropy of the electron distribution 
function in a collisionless plasma. We have adopted a two-dimensional velocity space Vlasov code in a simplified 
spatially one-dimensional configuration. 

We have found that the instability causes a violent deformation of the electron distribution function in phase 
space leading to the generation of short wavelength Langmuir modes. This mechanism coexists with the "fluid- 
like" generation of longer wavelength Langmuir modes due to the nonlinear modulation of the electron density. It 
corresponds to a non-local energy transport between modes in wave-number space. Eventually these short wavelength 
Langmuir modes lead to the formation of highly localized electrostatic structures. These structures correspond to 
potential jumps and are of interest for space observations. 
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